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Abstract 

We prove the existence of fractal solutions to a class of linear ordinary differential equations. This 
reveals the possibility of chaos in the very short time limit of the evolution even of a linear one 
dimensional dynamical system. 
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1 Introduction 



Understanding the complexities in a nonlinear dynamical system is of great interest in contemporary 
sciences. Both chaos and fractals seem to offer two important ingredients towards this effort. Even a 
simple quadratic nonlincarity as represented by the logistic map is known to generate, for a sufficiently 
large control parameter , a fractal attractor, indicating the onset of the deterministic chaos in the model. 
A linear equation, on the other hand ,is supposed to be immune to fractal attractors (solutions). In 
this paper we, however, report, on the contrary, a rather surprising fact: the simplest linear ODE of the 
form dy/dt = h{t)y admits fractal solutions, besides the standard exponential solution. This observation 
follows directly from our recent results indicating a deep relationship between the quantal geometric 
phase and the concept of time derived intrinsically from the (quantal) evolution [1] (denoted I henceforth). 
We give a brief review of this relationship in Sec. 2. In Sec. 3, we present the existence proof of the fractal 
solution for a (time dependent) first order linear ODE . 



2 Geometric phase and Time 

In a time dependent quantal evolution the evolving state acquires a geometric phase apart from the usual 
dynamical phase. Let the evolution of the quantal state 7p(t) be described by the Schrodinger equation 

iUj^ij = Hi; (1) 

where _ff is a time dependent Hamiltonian operator and t is the external Newtonian time. As time flows, 
the quantal state (< VlV" >= 1) traces a trajectory in the Hilbert space H, considered as a U(l) 
principal bundle over the projective space V. Due to the irreducible quantal uncertainty the actual path 
must always be fluctuating in nature thus allowing the evolving state to follow a path intersecting 'near 
by' rays in V. If the state is assumed to move along a cyclic path in V it returns eventually to the initial 
ray, but with a phase diff'erence. In the case of an open path, the phase difii'erence between the initial 
and final states equals to that of the cyclic path obtained by joining the initial and the final rays in V 
by the shortest geodesic. The total phase 7 so acquired by the state now consists of two components: 
7 — "fd + 7g, the dominant 'dynamical' and the relatively small 'geometric' phases respectively. The 
dynamical phase 7^ = J h{t)dt where h{t) =< il;\Hip > denotes the mean energy in the state, is a 
consequence of the mean dynamical evolution in the state. This mean evolution could further be realized 
as a pure verticle displacement of the state along the the fibre of its ray in the associated Hilbert bundle. 
Further, the scale of the external time t is set by the mean energy : t ^ -^,ho = h{to) , to=the initial 
time. As shown in I,the geometric phase -jg = J A where A = —i < (f)\d\(f) > and </> = e'^'^ip; on the 
otherhand could be interpreted as an effect of the inherent fluctuations in the actual quantal evolution. 
Although, the nature of this phase is of course geometric, as encoded in the associated parallel transport 
law, it appears ,in any way, as a correction term in the mean energy, thus offering itself to an equivalent 
dynamical (Hamiltonian) treatment. The geometric phase thus provides naturally a new (microscopic) 
time scale in the evolution, having a conjugate relationship instead with the energy uncertainty, rather 
than with the mean energy. Following the Leibnizian view [2], we call this time scale (variable) the intrinsic 
(geometric) time. Interestingly, this dynamical treatment of the geometric phase then allows one to write 
down a scaled Schrodinger equation which is suitable for an independent dynamical description of the 
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fluctuating state. The reproduction of a(n) (almost) self-similar replica of the Schrodinger equation in 
the scale of a fluctuation involves a sort of a 'renormalization' in the original equation which amounts 
to a subtraction of the mean (dynamical) motion in the state. It also turns out (see below) that the 
intrinsic time scale relates inversely (dually) to the extrinsic one in the limit of large fluctuations 0(1) 
(strong interaction) in the state. Thanks to the irreduciblity of the quantal uncertainty, one can then 
iterate the above steps an unlimited number of times, thus arriving at an interesting observation: quantal 
fluctuations are essentially fractals: having self-similar structures at all time scales. As a corollary, the 
trajectory of the evolving nonstationary state and hence the complex state function itself must also be 
fractals(at least in the time variable). The later conclusion was not stated explicitly in I. We shall supply 
a proof in favour of this in the next section. 

A remark is in order here. Quantum mechanics is a linear theory, with only a flrst order time 
derivative in the Scrodinger equation. Further, the (quantal) uncertainty principle is a direct outcome 
of this inherent linearity. Any linear time invariant model e.g., one in a signal processing system, also 
exhibits time-energy uncertainty. The above analysis should therefore be applicable not only in the 
domains stated above but also to any first order linear ODE with nontrivial time dependence. We recall 
that the (nonadiabatic) geometric phase was historically first discovered not in quantum mechanics but 
in optical polarisations having a linear Schrodinger equation-like guiding equation. 

3 Main Results 

Let us consider a first order linear ordinary differential equation of the form 

*^2/ = M% (2) 

where h{t) is a real function of the real parameter t. The factor i in the l.h.s. of eq(2) is introduced to 
keep the analogy with the Schrodinger equation in view. This is not necessary,however, in general. The 
complex functions y{t) are assumed to belong to a linear (Hilbert) space. The corresponding projective 
space thus consists of the rays y' = e*"y, a real. For the sake of clarity, we assume eq(2) to represent a 
linear time invariant dynamical system with mean 'energy' ho = h{to) at the intial time t = to- Eq(2) 
is thus (form) invariant under the group of translations, which in turn gaurantees the uniqueness of the 
standard exponential solution y = exp{—i J^* h{t)dt). We however wish to present a new class of fractal 
solutions to eq(2) by extending the group of translations to the (modular) group SL(2,R). We remark that 
the existence of a nontrivial geometric phase naturally provides the necessary window for introducing the 
duality transformation in a time dependent (linear) dynamical system. However, it turns out that the 
technique developed in I on the basis of a geometric phase is more general and could in fact be used in 
a situation with a vanishing geometric phase. 

Note that in a quantum mechanical model h{t) stands for a time dependent Hamiltonian operator. In 
an ordinary classical situation eq(2) may either be considered as a matrix equation e.g., the Fermi- Walker 
transport equation for the polarised light or be assumed to involve a background source of a motive force, 
to generate a nontrivial geometric phase. We however restrict the discussion to the simplest case where 
y denotes a single complex function having continuous first derivative. We show that SL(2,R) acts as an 
invariance group even in this apparently trivial equation. 

Let us proceed mimicking the steps of I thus indicating briefiy how an intrinsic time variable r appears 
in eq(2). Because of the time translation invariance, it is sufficient to study the nature of the solution 
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close to the initial time which we now set to t = 0. Further, assume f to be a dimensionless variable 
measured in the unit of ^.Thus the Taylor expansion of h{t) close to in eq(2) yields 

^^2/=[l + Ki)]2/ (3) 

The time dependent term in the rhs of this equation could be interpreted as the correction due to 'fluc- 
tuations' over the mean 'energy'. For a very short time interval (i.e.,neglictng O(f^) terms) , fluctuations 
can be approximated as u{t) = uit. In the external time scale fluctuations thus scale as t. However, v 

gives rise to a new time scale at the level of the first-order fluctuation, which follows from the equivalence 
of the description of the present model with that in I. Let us denote this intrinsic time scale by t ~ ^. 
To make the analogy with I complete let us write the 'Born-Oppenheimer' anstaz y = yoyi,yo = e~**, 
and introduce the intrinsic time r through 

di~ ^ di~ ~dT ^ ^ 

We thus obtain 

=vi[l + i^2(t)] 2/1 (5) 

Eq(5) has the following interpretation. The operator in the Ihs of (4) subtracts out the mean (dynam- 
ical) evolution in (3) that manifests in the unit of the external time 1). The residual (renormalized) 
evolution now consists only of the small scale fluctuations. However, in the scale of the intrisic time 
T ~ an externally small fluctuation does appear substantially large. Eq(5) then tracks this 'inter- 

nally large' fluctuation exactly in the same spirit as that of the original eq(3): i.e.; by identifying vi as the 
(scaled) mean 'energy' with a relatively small correction i'2{t) from the (2nd order) fluctuations. Further, 
the switching of the treatment from the external time scale down to the intrinsic scale amounts,in fact, 
to a transformation of the ordinary time (t) scale to a logarithmic scale: 

r = ln\t\-'^ (6) 

This is made explicit in the second equality in eq(4). Eq(5) thus offers a new approach in probing 
the very short time evolution of eq(2). Note that the standard exponential solution is obtained from the 
initial 'mean' solution yo near t = through a combination of two operations: i)by improving upon the 
initial solution with the inclusion of the neglected terms in the Taylor (perturbation) series expansion 
(the Picard's iteration), and then ii) by the successive applications of the translation group 'horizontally' 
along the ordinary t-axis. The first operation could indeed be achieved through the subtraction proceedure 
(c.f.,the Ihs operator in eq(4)), but continuing to work instead in the original time (scale) t. The new 
possibility that emerges in the present discussion is the foUowing.Instead of following the evolution further 
'horizontally', one could choose to dive 'vertically' down to the (intrinsic) logarithmic scale, thus climbing 
,as if, to a ripple of the background fluctuation. Consequently, the scale of the fluctuations gets sufficiently 
magnified (stretched), thus making a room for an (almost) independent treatment of the same, analogous 
to the mean motion. 

To continue, we note that the (-) sign in eqs(4) and (5) is typically a consequence of the geometric 
origin of the intrinsic time introduced in I. It reveals a sort of 'relativity' between the extrinsic and 
the intrinsic treatments: the direction of traversal of a path as seen from the external time frame gets 
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reversed in the intrinsic scale. Note also that the subtraction operation of the mean energy in every unit 
of time (e.g., f ~ 1) amounts to a folding( squzeeing) on the solution space. Eqs(3)-(5) thus encode a set 
of successive stretching -twisting and folding operations; which can clearly be iterated ad infinitum; thus 
establishing, in turn, a hierarchy of self-similar structures in the very short time limit of the evolution. 
A number of observations can now be made. 

1. Eq(6) represents a generalized duality relation between t and r. In the limit of small t, this reduces 
to i « t'~^,t' = 1 + r. Thus in the limit of large fluctuations, the SL(2,R) group is realized as the exact 
invariancc group of the ODE (2). 

2. The duality eq(6) also tells us that an exact replica of the original eq(3) is reproduced after the 
second iteration when t ^ 0. Succintly, the complete iteration process can be expressed as 

tn+l = -ln{Vntn), h = 1,1^0 = I, n = 0,1,2, .. . 

In the limit tn — > 0, the duality transformations reduce to the scaling relations tn+i = i^ntn- For, tn 
in the log scale tn+i = ^n{i^n^ — tn) and then replacing — tn by t„ (because of the translation 
invariance) one obtains the result. Thus the set of nonlinear operations, as detailed above, defines a 
hyperbolic iterated function system (IFS)[3] with scaling (contractivity) factor ui{= maa; (i^i, 1/2)), < 
ui < 1, at each point t = to oi the real t-axis. The intended fractal solution of eq(2) { &t t = to), our 
main result, is thus obtained as the imiquc attractor of this IPS. 

3. Interestingly, the limiting value of the scaling factor Vn, for any n, at the fixed point t„ = is given 
by the 'universal' value: Vg = , the golden mean. This follows from the necessary constraint that 
the system of equations (7) must coincide at the fixed point. One thus obtains the relations Vn+i{^ + 
Vn+i) = Vn^{l + Vn+2) = {i^n^ Vnli)i.K+i^n+2) ■ • ■ = 1> by virtue of the duality t„ i/^^ but Un+i ~ i„ 

4. As a consequencec,the limiting form of the set of iterated equations (7), at the fixed point, is given 

by 

i{d/dtn)yn = {-'^)"'Vgyn .gx 
tl = t, tn+l = Vgtn, n>0 

These equations could be interpreted as one obtained by splitting the infinite degenarecies of the original 
equation = y, through the repeated applications of the nonlinear stretching and twisted -folding 
transformations at the fixed point. Indeed, this could be realized by introducing a partition of the unity: 
1 = (1 — i>n) + Vn, identifying the braketed term as the 'mean' and the remaider as corrction due to the 
fiuctuations and then following the above steps. The apparent arbitrariness in the partition gets washed 
away at the fixed point, thus leading to the unique limiting system ,eq(8).0ne thus obtains eq(8) as the 
attractor of the class of ODE considered here, under the nonlinear invariance group SL(2,R). An explicit 
form of the fixed point fractal solution could therefore be constructed in the form 

y = exp{-i [t - Vg{ti -t2 + ta-.. .)]) (9) 

where the parameters t, tn are treated as independent variables and each tends to satisfying t>t\> 

t2 > t^ > By duality the point t„ = is mapped to tn+i = O,for each n, ensuring a nontrivial 

sewing of the ordinary exponentials along the (internal) 'verticle' direction. However,the function (9) 
collapses to the simple form y ~ exp{—iVgt) in the very short time limit, provided one makes use of the 
scaling relations tx = t,tn+i = Vgtn- 
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To interprete this result, we note that the original equation dy/dt = y,i = —it, in the log scale 
\i\ = lnT~^, can be translated as the definition of the box-counting dimension of the real axis parametrized 
by i : lim = 1, provided y is identified as the total number of infinitesimal intervals needed to cover 

T — ^0 

a finite segment of the f-axis.. In the present case, writing Ug = l — d,d= , we get instead a nontrivial 
scaling law y « T~^^~'^h Following I we interprete this result as one revealing a fractal structure in the 
real t axis itself. The real t axis thus behaves as a fat fractal with exterior dimension d[4]. The golden 

mean is thus realized as the corresponding uncertainty exponent of this fat fractal[l,4] . 

4. The above discussions also clearly reveal the presence of the deterministic chaos in the very short 
time evolution of a linear dynamical system. 

4 Final Remarks 

We have presented a new approach in analyzing a (linear) differential equation. The results discussed 
here are expected to find interesting applications in a number of physical problems where a perturbative 
approach normally fails. Ncvcrtlicloss, it's striking how the golden mean emerges as a universal scaling 
constant in probing the short distance (time) structure of a linear dynamical system. 
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